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Abstract. We demonstrate that by considering disordered single-particle 
Hamiltonians (or their random matrix versions) on ultrametric spaces one can 
generate an interesting class of models exhibiting Anderson metal-insulator 
transition. We use the weak disorder virial expansion to determine the critical 
value of the parameters and to calculate the values of the multifractal exponents 
for inverse participation ratios. Direct numerical simulations agree favourably 
with the analytical predictions. 
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A metric space is called ultrametric, if a distance function d(x,y) satisfies the 
strong triangle inequality d(x, y) < max{d(i, z), d(z, y)} for any three points x, y and 
z, so that all triangles are either isosceles or equilateral. The notion of ultrametricity 
appears naturally in p-adic analysis pQ which plays important role in various branches 
of mathematics and string theory. Ultrametric spaces have found many useful 
applications in physics and beyond [2], most notably in the theory of spin glasses [3] 
where they describe the hierarchical organization of free energy landscapes. Random 
matrices whose structure reflects some specific aspects of ultrametricity were under 
active study recently due to their relevance to description of relaxation properties in 
complex landscapes typical for heteropolymers and related systems [4]. 

The aim of the present work is to demonstrate that by considering disordered 
single-particle Hamiltonians (or their random matrix versions) on ultrametric spaces 
one can generate an interesting class of models exhibiting Anderson metal-insulator 
transition. The latter is one of the central phenomena in the theory of disordered 
quantum systems and describes an abrupt change in the nature of eigenstates from 
extended to spatially localized. Exactly at the transition point eigenstates are very 
inhomogeneous and strongly fluctuating and can be characterized by non-trivial 
multifractal dimensions. Despite a great progress in understanding of the Anderson 
transition achieved in the last decades |5] based mainly on the mapping to the nonlinear 
it— model [5J [7], there are still relatively few models of disordered Hamiltonians 
exhibiting the localisation transition and affordable for direct analytical studies. 
Among those deserving special mentioning are the models on the infinite tree-like 
graph (Bethe lattice) [8] and the ensemble of power-law banded random matrices 
introduced in [9] and further investigated in [lOj which became rather popular in recent 
years. The problem of rigorous mathematical analysis of the Anderson transition is 
yet outstanding, see e.g. [TT] . 

The idea that models of statistical mechanics on hierarchical spaces with 
effectively ultrametric structure can be efficiently studied analytically due to their 
high symmetry goes back to Dyson 12J and proved to be rather useful in shedding 
light on renormalisation group methods and developing rigorous mathematical proofs 
in general theory of phase transitions |13) . Recently there was a systematic attempt 
to put disordered models of spin-glass type on such lattices [14]. We hope that the 
models proposed in the present paper will provide a useful framework for further 
insights into the nature of the Anderson localisation transition, and perhaps could 
yield to a rigorous mathematical analysis. 

Let us start by introducing a random matrix ensemble of Hermitian matrices 
with ultrametric structure whose eigenvectors exhibit the Anderson transition when 
changing some control parameter. The idea behind such an approach is that ensembles 
of that sort can be not only efficiently studied numerically, but should be also amenable 
to a variety of existing analytical techniques like, for example, mapping to a kind 
of nonlinear <r— model on ultrametric space, which could keep promise for an exact 
solution (see a short discussion in the end of the paper). Leaving the investigation of 
these issues to a future research, in this paper we are to exploit some recent analytical 
insights provided by works on almost-diagonal random matrices [HI [16] . 

The ultrametric structure for a particular representative H of the suggested 
ensemble of random matrices is encoded in the choice of variances for the Gaussian- 
distributed entries. Namely, all 2 K x 2 K entries Hij of the matrix H are taken 
as independent complex Gaussian random variables (up to the obvious Hermiticity 
constraint H*, = Hji) with the zero mean value. All diagonal elements of H are 



Anderson localization transition in ensemble of ultrametric random matrices 



3 



(a) 




Figure 1. (a) Schematic structure of variances of entries of an ultrametric random 
matrix, (b) Hopping between nodes of the ultrametric lattice. 



real and have the same variance: (H^) — W 2 , Vi. The variances of the off-diagonal 
elements have a block structure typical for " hierarchic" lattices (Fig. QJ,) , with blocks 
of the sizes 1, 2, 2 2 , . . . , 2 K ~ 1 . Inside each block of the size 2 l the variances (iJ^i/ij) 
are constant and chosen to be J 2 /p 21 , controlled by two real parameters J and p. 
Such a random matrix describes random hopping between boundary nodes of a tree 
of K generations with given coordination number which we choose for simplicity to 
be equal to two (Fig. [TJd) . The underlying tree serves the purpose of defining the 
ultrametric distance d(i,j) between any pair of boundary nodes. It is defined as the 
number of edges in the shortest path connecting nodes i and j. As easy to check 
the corresponding distance is ultrametric. The formula for the variances can be then 
expressed naturally in terms of d(i,j) as (H^Hij) = J 2 /p d ( l >^~ 2 . 

The ultrametric random matrix model introduced above exhibits the Anderson 
transition by changing parameter p. Namely, in the thermodynamic limit K — > oo the 
eigenvectors of H are extended for p < 2, localized for p > 2 and critical (multifractal) 
at p = 2 . We expect this scenario to hold for any ratio J/W, implying in particular 
that at p — 2 changing J/W should allow to pass from strongly to weakly multifractal 
eigenstates. Assuming that J/W <C 1, the fractal dimensions d q of the eigenvectors 
(defined below in Eqs.([2][3|)) in the critical regime can be calculated explicitly and are 
given by: 



J r( 9 -i/2) 



d a = 
q WV21n2 



q > 



(i) 



T(q) ' " 2' 
where T(q) is the Euler Gamma- function. 

We define the moments of the eigenvector components I q (also known as inverse 
participation ratios, IPR's) in the standard way: 

1 N 

n=l 



P q (E,n) 



P(E) 



N 



* k {n)\ 2 «5(E-E k ) 



(2) 



where the normalized eigenvectors ^ with ./V = 2 K components ^.(n), n = 
1,2, ... ,N correspond to the eigenvalues E k and p(E) stands for the mean density of 
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states at the energy E (which will be fixed to zero for simplicity). At criticality the 
ensemble-averaged moments of the eigenvectors scale anomalously with the system 
size N, defining a set of fractal dimensions d q : 

(I q ) ocN-^-V. (3) 

In order to calculate the moments of the eigenvectors under the condition 
J/W < lwe use the method of the virial expansion developed in [16l[T7]. The methods 
allows us to construct a well-controlled perturbation expansion of the moments of the 
eigenvectors using the small parameter J/W. The first term of the expansion is 
generated by the diagonal part of a matrix H . All corresponding eigenvectors have 
only one non-zero component, which is equal to one due to the normalization condition. 
As a result all moments of the eigenvectors are equal to one in this approximation: 
'q°^ = 1. The corrections to this result originating from the off-diagonal matrix 

elements of H can be represented as certain integrals over commuting and anti- 
commuting variables. That integrals can be either calculated explicitly [16[ I17j or 
alternatively they can be interpreted as the moments of the eigenvectors of certain 
submatrices of a matrix H. To be specific, the lowest order in J/W correction Iq is 
generated by the following 2x2 sub-matrices of H : 

tj t \ ( Hnn Hnm \ -i at / 

n[n,m)= „„ „ I , n, m = 1, . . . , JS, n^=m 

\ "tim H-mm ) 

Denoting by P q (n,m) the moments of the first component of the corresponding 
eigenvectors we are able to represent (iq 1 ^ as follows: 

1 N 

I P) = x y E« p <,( n > m ))- 1 ) ( 4 ) 

The subtraction of unity is required in order to get rid of the contribution already 
taken into account at the previous step. The above approximation implies that the 
pair- wise "interaction" of the levels of the unperturbed system is taken into account 
at this stage. Essentially the same idea was underlying the real-space renormalization 
approach developed for critical systems with a long-range interaction and the power- 
law random banded matrix model 1181 1101 . 



Below we calculate the averaged value of I q using an explicit expression for 
the eigenvectors components. To this end let us introduce the following notation 
Ei = H nn , E% = H mm , h = H nm and P q = P q (n,m). The eigenvalues Ai, A2 of the 
matrix H (n, m) are equal to 

1 



A ll2 = -(E 1 + E 2 ± - E 2 f + 4\h\ 2 ). (5) 

The absolute values of the components ^1.2 of the corresponding normalised 
eigenvectors are given by libi ol = = Averaging the moments of \ibi\ 

6 6 J \Y1,1\ y/\h\ 2 + (\-E 1:2 ) 2 6 6 IV 1 

first over Ex and E 2 we obtain the following expression for (P q ) 



E- 



{ q ' E ~ J ./-co 2*PW 2 (\h\? + (A; - E,fY ' W 

where p is the density of states at E = 0. The latter can be easily calculated in the 
diagonal approximation giving the result p = l/y/2nW. The corrections generated by 
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off-diagonal elements of H is of the order of (J/W) 2 [TS] and hence can be neglected. 
Eq.® shows that A x = if E X E 2 = \h\ 2 and E 1 +E 2 < 0, while A 2 = if E X E 2 = \h\ 2 
and Ei + E 2 > 0. Therefore we should consider separately two domains of integration 
E\ + E 2 < and E\ + E 2 > 0. In each domain only one of the 5-functions in Eq.© 
survives. After a simple change of variables we are able to integrate over E\ in Eq.([6]): 

l''l 4 + E l / rp2 \ I- 1 



P n „ = / ' e 



q ' E J-ocV^W ' \E 2 + \h\ 2 



J 2 



dx\h\ _\h^^ +x ^ 



2 \ 9- 1 



e-***^ * -f— , (7) 



<2^W \x 2 + l / 

where the new variable of integration x = \h\E 2 is introduced. Now we are ready to 
perform averaging over the complex variable h. Since the expression above depends 
only on \h\ the integration over the argument h is trivial and just leads to appearance 
of the factor 2ir while integration over \h\ can be done explicitly: 

(p) — — r dx (8) 

z J-oo (1 + a (x A + x 1 )) 1 

where a = \/{\h\ 2 )/ \2W . The asymptotic behaviour of this expression in the limit 
a — > has the following form: 

- r(g-i/2) 
r(?-i) 

The first term in the expansion corresponds to the diagonal approximation and is 
cancelled by —1 in Eq.((3]), while the second one gives a non-trivial contribution. 
Collecting the results of Eqs. (|2l4l9|) we find (I q ) up to the first order in J/W: 

N 



(P«) E ,h = 1 - ^ W a + °(« 2 )- ( 9 ) 



/F r( g -i/2) i A 
w 1 V2 r( 9 -i) ]v4- ^ (10) 

= J /F r(g-l/2) " 1 fnl 

T^V 2 r(g-l) ^ pd(l,m)/2-i • v ; 

According to the definition of the ultrametric distance d(n, m) the sum in Eq. is 
given by the geometric series: 

/ \ K 

A_ i ^/o\' (I) -1 



s - Z p d(i,m)/2-i - Z ( p ) - v a - 1 ■ ( 12 ) 

m=2 ^ i=0 V ^ 7 P 

As if — > oo the asymptotic behaviour of the above sum crucially depends on the value 
of the parameter p: for p > 2 the sum tends to the finite value S — > p/{p — 2) < oo 
independent of A implying that eigenvectors are localized, whereas in opposite 
situation p < 2 we see that S diverges as S « (p/(2 — p))A 1_log2 p , signalling of 
a breakdown of the perturbation theory for arbitrarily weak off-diagonal coupling. 
It is natural to interpret such breakdown as the signal of eigenvector derealization. 
Finally, at the critical point p = 2 we have S = K = log 2 N leading to 

Comparing this result with the definition of d q given by Eq.Q, and expecting 

for J/W <C 1 we interpret this marginal behaviour in terms of the anomalous scaling 
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of the moments of the eigenvectors, and the expression for the fractal dimensions (TTJ) 
immediately follows. The condition q > 1/2 is necessary to guarantee the convergence 
of the derivative with respect to a of the integral in © at a — 0. One expects 
that the region q < 1/2 can be covered by employing the powerful symmetry relation 
discovered in [19] , see below. 

One can notice that the g-dependence of the above expression is exactly the same 
as the one found for the ensemble of the power-law random banded matrices [10 . The 
matrix elements of that ensemble are independently distributed Gaussian variables 
with zero mean and with variance decaying in a power-law fashion as a function of the 
distance to the main diagonal. Such a coincidence can be easily explained by Eci. tfTO"]) . 
Indeed, we can see that g-dependence enters only in the pre-factor and thus is always 
the same for any Gaussian ensemble with independent entries. However we would like 
to stress that this " superuniversal" behaviour should be expected only in the regime 
J/W < 1. 

The results of numerical simulations confirming our predictions are presented in 
Fig. [2] The moments of the eigenvectors, whose eigenvalues are close to E = 0, were 
calculated for various system sizes ranging from N — 2 s to N = 2 14 . To obtain 
their averaged values reliably a large number of random matrices were generated, so 
that the number of eigenvectors contributing to the computed averaged values were 
ranged from 10 3 for N — 2 14 up to 10 5 for N = 2 8 . The extracted values of the fractal 
dimensions d q for J = 1, W = 10 and W — 1/3 are represented in Fig. [2] by circles and 
triangles respectively. For q < the eigenstates intensities were first coarse-grained 
over boxes containing 8 sites. The values of d q for W = 10 are in perfect agreement 
with our analytical prediction (solid line in Fig. [2]). The latter was calculated directly 
from Eq.([T]) for q > 1/2. To make the prediction for q < 1/2 we used the remarkable 
symmetry relation discovered in [19] and intimately related to profound fluctuation 
relations in the field of non-equilibrium dynamics 20J . According to it the anomalous 
multifractal exponent defined as A q = (d q — d)(q — 1), where d is the dimensionality 
of the space, is symmetric with respect to line q = 1/2, i.e. A q — Ai_ g . In our model 
d = 1, as according to Eq.© the fractal dimensions are defined by the scaling of the 
moments of the eigenvectors with respect to the total number of sites N ("volume" 
of the system). We checked the validity of the relation for our model (inset of Fig. [2]) 
and employed it to calculate d q for q < 1/2 analytically using Eq.([T]). The fact that 
the expression (p} diverges at q = 1/2 indicates that our perturbative approach breaks 
down at some values close to q = 1/2, at which Eq.fl]) becomes of the order of one. 
However we would like to mention that the values q* can be made arbitrary close to 
1/2 by taking J/W sufficiently small. 

Another possible realization of the ultrametricity in the framework of disordered 
systems is the ultrametric Anderson model, which is a version of the Anderson model 
defined on the lattice of N = 2 K sites supplied with the ultrametric distance d(i,j) 
described above. The diagonal matrix elements of the Hamiltonian are assumed to 
be random and characterized by the variance W. The off-diagonal matrix elements 
Hij of the Anderson model are however non-random. Their values = J/p d ( l >-?)/ 2_1 
depend on two parameters J and p and the distance d(i,j) similar to corresponding 
variances of an ultrametric random matrix. 

Although one can not apply directly our results to the Anderson model, one 
can use the argument due to Levitov [T5] to establish the existence of the Anderson 
transition for this model. According to that argument, two sites i and j will be mixed 
by the off-diagonal entries of TC only when they are " in resonance" : \Hu — Hjj | < . 
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Figure 2. Fractal dimensions d q for J = 1, W = 10 (circles) and W = 1/3 
(triangles). The solid line is the analytical prediction based on Eq. JTJ. Inset: 
numerically calculated exponents A q (solid) and Ai_ 9 (dashed). The gray shaded 
background marks the 95% confident interval of the numerical results. 



It is easy to estimate the probability that two sites are in resonance: 

= / / ,!! ^±< ' " ' JU - \Hu H n \) (14) 

where 6(x) is the step-function. Calculating the integral yields P« = 
jH tj /V2W e _ x /2 ^ T ^ ^/l"^. m our approximation. Then the mean total 
number of sites which are in resonance with a given site i is given by: 

When the number is divergent in the limit N — > oo as a power of iV the couplings, 
whatever weak, produce derealization of eigenvectors in the site basis. On the other 
hand that number tending to a Unite limit indicates localization of the corresponding 
eigenvectors in the weak coupling limit (which of course does not preclude possibility 
of derealization at some critical coupling). Finally, the transition between the two 
regimes happens precisely when {4f) i ~ In AT. For the ultrametric Anderson model 
T,j 7 a H ij/W = {J/W)S, where S is defined by Eq.JTJ]). Thus we conclude that 
depending on p the nature of the eigenvectors of the Anderson model is the same 
as for its random matrix analogue: localized for p > 2, extended for p < 2 and 
critical at p = 2, although we expect the dependence of anomalous dimensions d q 
on the parameter J/W to be different in the two models. The statistics of energy 
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level spacings is expected to be Poisson for p > 2 and Wigner-Dyson for p < 2. For 
earlier studies of spectral properties of the random Schroedinger-like operators with 
hierarchical structure see [STJ [55] . 

Finally, let us say a few words about possibilities to understand the opposite 
limit of strong off-diagonal coupling J 3> W. In the original formulation of the model 
neglecting diagonal entries does not seem to lead to any essential simplification, and 
we are unfortunately unable to provide any insights to the problem. However, one 
can think of introducing an "n-orbital" variant of our model in the spirit of [6]. 
In this way one should replace each entry Hij of the hierarchical random matrix 
Hamiltonian with n x n block of independent, identically distributed zero mean 
Gaussian elements with the variances J 2 /p 21 inherited from the original n = 1 model. 
Letting n — ► oo and simultaneously rescaling the parameters J and W in such a 
way that J 2 /W 2 = t/n with some t — 0(1) ensures a non-trivial limit of the theory 
equivalent to a lattice version of the standard nonlinear er— model, with couplings 
depending on the ultrametric distance between the lattice sites. In such approach the 
strong coupling limit t 1 will be amenable to the standard perturbative expansion in 
diffusion- like modes like in [9] , with the diffusion propagator replaced by its analogue 
on the hierarchical lattice (note that as the model lacks the translation invariance the 
standard Fourier modes expansion is ineffective and should be replaced by a certain 
ultrametric analogue). We leave analysis of this and related questions to subsequent 
publications. 
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